Phase-slippage and self-trapping in a self-induced bosonic Josephson junction 
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A dipolar condensate confined in a toroidal trap constitutes a self-induced Josepshon junction 
when the dipoles are oriented perpendicularly to the trap symmetry axis and the s-wave scattering 
length is small enough. The ring-shaped double-well potential coming from the anisotropic character 
of the mean-field dipolar interaction is robust enough to sustain self-trapping dynamics, which takes 
place when the initial population imbalance between the two wells is large. We show that in this 
system the self-trapping regime is directly related to a vortex-induced phase-slip dynamics. A 
vortex and antivortex are spontaneously nucleated in the low density regions, before a minimum 
of the population imbalance is reached, and then cross the toroidal section in opposite directions 
through the junctions. This vortex dynamics yields a phase slip between the two weakly linked 
condensates causing an inversion of the particle flux. 

PACS numbers: 03.75.Lm, 03.75.Hh, 03.75.Kk 



Phase-slippage is a dynamical process that generally 
takes place in ac Josephson junctions. In this regime the 
phase difference increases linearly in time and propor- 
tionally to the external voltage applied or chemical po- 
tential difference, respectively in superdonductor or su- 
pcrffuid junctions. This linear growth can be interpreted 
as a periodic change of the phase difference by an amount 
of 27r, or as the system periodically undergoing a phase- 
slip. The current understanding is that phase slips are 
related to the dynamical creation of vortices, which cross 
the flow path and leave a 2ir phase behind PJ ■ This pro- 
cess has been widely addressed in supcrfluid helium (see, 
for instance, 0) and is receiving increasing attention in 
Bose-Einstein condensates (BECs) 0-01 ■ 

In the field of supcrfluid He it has been shown that 
after the phase-slip takes place the superfluid velocity is 
reduced by a quantized amount, and under some condi- 
tions the flow of atoms can even be reversed. The condi- 
tions for reaching the phase-slip regime depend strongly 
on the geometry of the constriction, that is the weak link 
that creates the Josephson junction. When the width of 
the junction is larger than the healing length then phase 
slippage sets in and the current-phase relation is linear, in 
contrast to its sinusoidal shape in the Josephson regime. 
The transition from phase slips to Josephson regime has 
been observed in superfluid 4 He by increasing the healing 
length with respect to the size of the constriction [(| . 

In BECs, Josepshon junctions are experimentally 
achieved by confining the atoms in an external double- 
well potential. In these systems the observation of 
plasma and self-trapping oscillations has been reported 
[3, as well as the ac and dc Josephson effects Q. Also, 
Josephson oscillations in an array of junctions have been 
addressed |{|, demonstrating phase coherence between 
them. On the other hand, phase-slips in BECs have been 



recently investigated both theoretically [H, H| and experi- 
mentally Q as a dissipation mechanism of the superfluid 
flow in toroidally confined condensates with a constric- 
tion produced by an external laser beam. 

Constrictions to the superfluid flow need not be cre- 
ated by an external potential, but can also be brought 
about by the condensate itself. In particular, in a dipolar 
Bosc-Einstcn condensate (dBEC) confined in a toroidal 
trap, the anisotropic character of the dipolar interaction 
induces two potential barriers along the torus, which can 
be regarded as constrictions [loj . They can be achieved 
by polarizing the dipoles perpendicularly to the trap sym- 
metry axis and tuning the scattering length to small val- 
ues. Under these conditions, a self-induced double- well 
effective potential appears in the system, due to the com- 
bination of the toroidal trapping potential and the mean- 
field dipolar interaction. In the self-induced double-well, 
a dBEC behaves as a self-induced Josephson junction 
(SIJJ) [ll|, showing both Josephson and self-trapping 
(ST) regimes. In this letter we study in detail the self- 
trapping regime and show that it is closely related to 
phase-slip dynamics. 

In the mean-field framework, the time evolution of 
the system is given by the 3D time-dependent Gross- 
Pitaevskii (TDGP) equation: 
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where 'J = \I/(r,i) is the condensate wave function, m 
is the atomic mass, and g = 4iTh 2 a/m is the coupling 
constant of the contact interaction, with a the s-wave 
scattering length. The toroidal confinement is simu- 
lated by a harmonic plus Gaussian potential: Vt(r) = 
(to/2) w 2 (r^ + \ 2 z 2 ) + Vq exp(— 2r 2 _/ Og) , where ujj_ is 
the radial harmonic trap frequency, A = uj z /uj± is the 
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trap aspect ratio, with u> z the frequency in the z direc- 
tion, and Vo and do are the intensity and waist of the 
laser beam producing the central hole of the torus. For a 
condensate of polarized magnetic dipoles, the mean-field 
dipolar potential can be written as 
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where [Iq is the vacuum permeability, /i the magnetic mo- 
ment of the atoms, and 9 is the angle between r — r' and 
the magnetization axis. The wave function ^ is found as 
a function of time by evolving Eq.([T]) in real time. We ap- 
ply Hamming's algorithm (predictor, corrector, modifier) 
initialized with a fourth-order Runge-Kutta method, and 
treat the dipolar term using Fourier transform techniques 
(see [ll[ for further details). 

In the calculations we consider a condensate with 
N = 5 x 10 4 atoms of 52 Cr with a = 14 as and fi = 6/xb 
(hb is the Bohr magneton), confined in a toroidal poten- 
tial with wj_ = 8.4 x 2tt s _1 , A = 11, Vq = 30hcj^ and 
(To = 2.08 ax, where a± = ^Jh/muj^ is the radial har- 
monic oscillator length. Since the system is scalable, it 
can be map ped into other trap parameters and bosonic 
species [1JJ. We want to stress that the trap symmetry 
axis lies along the z axis, whereas the dipoles are oriented 
perpendicularly to it, along the y axis. In this configu- 
ration the effective potential V e s(r,t) — Vt(r) + Vd(r,i) 
is a ring-shaped double well with the two wells on the 
x axis. Note that this effective potential depends on 
time through the condensate wave function in the dipolar 
term. The self-induced double-well structure is not ex- 
clusive of dipolar condensates [l(| , but has also been pre- 
dicted to appear in degenerate dipolar Fermi gases [12j . 
and recently in more general dipolar Bosc and Fermi sys- 
tems [13l | . However, coherent tunneling dynamics is only 
possible if the system is in the supcrfluid regime. 

Dynamics in the SIJJ consists on the coherent tun- 
neling of atoms between the two weakly linked dBECs. 
It is generally well captured by two conjugate variables: 
the population imbalance, Z(t) = (N L (t) - N R (t))/N, 
and the phase difference, </>(t) = 4>n(t) — 4>L(t), where 
NL(R)(t) corresponds to the number of atoms on the left 
(right) well, and 4>l(r) (t) to the phase of the dipolar con- 
densate averaged on the left (right) well. Note that this 
description assumes, in particular, that the phase of the 
wave function is uniform in each well. For narrow junc- 
tions this is a good approximation, but it starts failing 
for wider junctions. 

For initial conditions Z(0) = 0.67 and 0(0) = the sys- 
tem is in the self-trapping regime. This regime is brought 
about by the nonlinear interactions in the condensate 14 1 
and is characterized by a nonzero time average of the pop- 
ulation imbalance, that is, the trapping of the most part 
of the atoms in one of the wells. It is reached for strong 
enough interactions and for a large initial population im- 
balance between the two wells. In the cases observed 
experimentally Q, ST was accompanied by a linear in- 
crease of the phase difference between the two weakly 



linked condensates, characterizing thus a running-phase 
mode. 

Figure [Ha) shows the effective potential at t = 0, in 
the z = plane. It has the shape of a double well in 
the azimuthal direction, with two potential barriers at 
y > and y < 0, which we will name as upper and 
lower junction, respectively. The two junctions are cou- 
pled in phase, so they can be effectively regarded as only 
one, analogously to the array of nanoapertures in helium 
experiments (see, for instance, @). The initial density 
structure with the two links is shown in Fig. [ljb). As a 
guide to the eye the dashed line draws the equidensity 
contour at 0.1/9 max , being p max the maximum density in 
the plane z = 0. We have taken into account in the 
definition of left and right wells the displacement of the 
junctions around x ~ a± with an inclination of about 
15° from the y axis. This gives rise to a left well slightly 
larger than the right one [15| . 
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FIG. 1: (Color online) Initial effective potential (a) and cor- 
responding density (b) in the plane z = 0. Imbalance and 
phase difference as a function of time. 

Fig. (He) shows the coupled dynamical evolution of 
Z(t) and 4>(t) for the above initial conditions. The time 
derivative of the population imbalance is related to the 
flow of atoms across the junction. At the maxima and 
minima of the imbalance its slope changes sign, so that 
the flow of atoms is reversed. In a two mode description, 
when the imbalance reaches a maximum the phase differ- 
ence is zero, while in a minimum it reaches n [l6j . How- 
ever, looking more closely to the coupled dynamics given 
by TDGP, we see that in the minima of imbalance the 
phase difference does not exactly correspond to 4> = n. 
This fact hints at a richer dynamics that is not accounted 
for in the two macroscopic variables Z and <j), but comes 
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from a local behavior. In other words, since the junction 
is a wide junction, the local currents within it affect the 
gross dynamics contained in Z and <j). To understand bet- 
ter the irregularities in the signal of Fig. [TJc), we have to 
look at the dynamics of the local phase. We find two main 
results. On the one hand, we see that the flux inversion in 
the minima of the imbalance is accompanied by the cross- 
ing of a vortex through each junction. These vortices are 
dynamically created and constitute the phase-slip pro- 
cess that takes place when the phase difference reaches 
7T. On the other hand, the currents generated by the 
phase gradients within the junction propagate and affect 
the phase coherence in each one of the wells, giving rise to 
a spatial dependence of the phase. This new phase gra- 
dients become currents in the wells that are not properly 
accounted for in Z and <j> and are mostly the responsible 
of the dephasing between these two variables. 

Let us examine the dynamical nucleation of vortices 
in more detail. In the top panels of Fig. [2] we schemati- 
cally describe the vortex-induced phase-slip process that 
contributes to the flux inversion in the minima of the im- 
balance. Before reaching the minimum, since the phase 
difference is very large, there is an accumulation of phase 
between the two condensates, producing a large gradient 
of the phase along the line that connects both junctions. 
The atoms tunnel from the left to the right well, acquir- 
ing a velocity associated to this gradient, as indicated in 
Fig.[5Ja). The velocity is larger in the low density regions, 
namely the central hole giving rise to the toroidal geom- 
etry and the external region. When the phase difference 
is 7r, the velocity reaches a critical value and two vortices 
are nucleated, a vortex-antivortex pair to ensure angular 
momentum conservation, as seen in Fig. [UJb) [l7| . Sub- 
sequently both vortices cross simultaneously the section 
of the torus in opposite directions, each one moving out- 
wards through one link, producing the phase slip in the 
junction. After the passage of the vortex (antivortex) 
the sign of the local velocity in the link is changed, lead- 
ing thus to the inversion of the flux of atoms, as shown 
schematically in Fig. [He). The vortex (antivortex) with 
positive (negative) charge is defined with the phase gra- 
dient in the counterclockwise (clockwise) direction. 

Panels (d), (e) and (f) of Fig. [2] show the phase of the 
wave function, in the z = plane of the upper junc- 
tion, around the first minimum of Z(t) at the moments 
schematically represented in (a), (b) and (c), respectively. 
A vortex-antivortex pair is nucleated in the center, sep- 
arates and then the vortex (antivortex) crosses the lower 
(upper) junction. The same effect is produced if an an- 
tivortex and a vortex are nucleated at opposite points 
of the external surface of the condensate and then cross 
the torus and anihilate at the center. Analyzing the dy- 
namics we see that both processes take place at different 
minima. A similar phenomenon has been addressed in 
Ref. Q, where a single external barrier was raised in a 
toroidally confined s-wave condensate with a nonzero ini- 
tial angular momentum. In that work the appearance of 
vortices from the external or internal surfaces of the torus 
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FIG. 2: (Color online) A scheme of the slippage process is 
given in the top panel. The solid lines represent constant 
phase curves, the dashed lines a torus contour in the z — 
plane. The + (— ) signs the position of a vortex (antivortex), 
and the arrows indicate the phase gradient, and thus the ve- 
locity field, along the junctions. We have drawn the following 
three steps: a) Phase accumulation along the junctions, b) 
Vortex-antivortex pair separation after the nucleation and c) 
Inversion of the velocity field after the vortices cross the torus. 
In the bottom panel we show snapshots of the numerical phase 
around the upper junction at z — 0, as a function of x and 
y. The times are: (d) t = 6.41 wj 1 , (e) t = 8.12uj 1 , and (f) 
t = 8.76 W7 1 . 



was controlled by the height of the barrier, while the 
nonzero angular momentum favored the first process. In 
the case reported in this article, however, both processes 
are equivalent. We want to note that all the vortices 
we predict correspond to quantized vortex lines, with no 
visible bending along the z direction. The quantization 
is guaranteed by the fact that the phase changes by 2ir 
around the singularity. The angular momentum remains 
zero throughout the simulation and there is no net cre- 
ation of vorticity in the system. 

In order to complete the description of the vortex- 
induced phase slip, we have computed the mean velocity 
of the atoms across the upper and lower junction as a 
function of time, see Fig. [3] As expected, both junctions 
behave in phase and present the same physics: simul- 
taneously to the nucleation and crossing of the vortices 
through the junctions, there is a steep drop of the veloc- 
ity in each junction. In analogy with the experiments 
carried out in supcrfluid helium, this drop marks the 
point at which a phase slip takes place. From the sharp 
decrease of the mean velocity it is possible to estimate 
the vortex passage velocity using v p = A./ St, where A 
is the length of the barrier and St is the time it takes 
the vortex to cross it as read from Fig. [3] For the three 
minima it gives: Vp = 2.7 u;±a±, v p 2 ^ = A.2oj±a± and 
= 8.2 u!±a±. 
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FIG. 3: (Color online) Velocity of the atoms at the upper 
and lower junctions as a function of time. The steep drops 
correspond to the crossing of a vortex (phase-slip). As a guide 
to the eye, the imbalance has been plotted in an arbitrary 
scale as a dashed line (see Fig[TJ . 

In conclusion we have shown that there is a close rela- 
tionship between phase-slips and self-trapping dynamics 
in a ring-shaped Joscphson junction. The coherent pas- 
sage of a vortex and an antivortex across the two coupled 
junctions triggers the inversion of the flux of atoms in the 
minima of imbalance, where the phase difference reaches 
7r. These results are the first step to bring together the 
processes of self-trapping, ac Josephson effect and phase- 
slippage and to gain a deeper understanding of non-linear 



dynamics in bosonic Josephson junctions. 

We believe that these physics can be experimentally 
addressed, since the needed techniques are nowadays 
available: toroidal traps are currently constructed in ac- 
tual experiments flU [l9| . The initial imbalance can be 
achieved using cither an external ramp or by shifting the 
laser beam. The reduction of the scattering length is 
provided by the utilization of Fcshbach resonances pro- 
cedures. Our results involve Cr atoms [20j], but also Li 
or K are candidates to be experimentally studied plj . 
Both density and phases have been measured in experi- 
ments Q ■ Experimental observation of vortex-dipoles in 
Bosc-Einstein condensates [22| has also been addressed. 

Finally we want to point out that, although we have 
studied the case of a SIJJ, the relation between self- 
trapping and vortex dynamics is not restricted to dipolar 
condensates but it can be extended to any other type of 
bosonic Josephson junctions [23| . 
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